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Abstract. - We study the optical loading of a trapped Bose-Einstein condensate by spon- 
taneous emission of atoms in excited electronic state in the Boson-Accumulation Regime. We 
generalize the previous simplified analysis of ref. [Phys. Rev. A 53, 2466 (1996)], to a 3D 
case in which more than one trap level of the excited state trap is considered. By solving the 
corresponding quantum many-body master equation, we demonstrate that also for this general 
situation the photon reabsorption can help to increase the condensate fraction. Such effect 
could be employed to realize a continuous atom laser, and to overcome condensate losses. 
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During the last years the experimental realization of Bose-Einstein condensation (BEC) 
Oh, in trapped weakly interacting gases 0] has stimulated an enormous interest [Q. Among the 
results related to BEC, specially remarkable is the realization of an Atom Laser As 
a coherent source of matter waves, the atom laser will lead to new applications in atom 
optics. Several groups have demonstrated atom lasers using pulsed H,|(| or continous || 
qh outcoupling from the BEC, either by employing rf fields ||, or Raman pulses ||. However, 
the continuous outcoupling is just a halfway towards a cw atom laser. The continuous loading 
of the condensate still remains to be realized. Like in the developement of light lasers the 
availability of cw atom lasers would open the way to "high power" and precision applications. 

Two different physical mechanisms could provide continuous pumping of atoms into a 
condensate. On one hand, this can be achieved by collisional mechanisms fij, in which two 
non-condensed atoms collide, one being pumped into the condensate, whereas the other being 
evaporated. On the other hand, the optical pumping into a BEC via spontaneous emission 
has also been proposed ||. If the reservoir could be filled in a (quasi-) continuous way by 
laser cooling techniques, one would benefit from the large cooling efficiency of laser cooling 
compared to evaporative cooling, allowing for a considerable increase in atomic flux produced 
by an atom laser. For the latter, it is crucial that the spontaneously emitted photons cannot 
be reabsorbed, since otherwise heating is introduced in the system, and BEC can be neither 
achieved, nor maintained It has been shown that in the so-called Festina Lente regime, in 
which the spontaneous emission rate 7 is smaller than the trap frequency lo, the heating due 
to the reabsorption processes is suppressed |Tc| ]. However, due to the slow time constants in 
this approach the cooling efficiency is greatly reduced. Recently, it has been also reported [JOl 
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that when an atom possesses an accessible three level A scheme, where one of the atomic 
transitions decays much faster than the other, the reabsorptions in the slow transition can be 
largely suppressed, without the time limitations of Festina Lente. 

In this Letter we study the loading of a BEC, which is formed in a trapped electronic 
ground-state \g), via spontaneous decay of atoms from an also trapped internal excited state 
\e). In the regime Nq 3> a, N — No (where a is the effective number of levels other than the 
condensed one to which the excited atoms may decay, N is the total number of atoms in the 
\g) trap, and 7Vo is the number of condensed ones) the atom will decay into the condensate 
with high probability, due to the bosonic enhancement. Such regime has been named Boson- 
Accumulation Regime (BAR) [ |l2| . Note, however, that even in the BAR, in the absense of 
reabsorption simple arguments imply that the mean proportion of atoms in the condensate 
after the decay decreases |^2[ . If the reabsorptions are present, each subsequent decay would 
lead to an even larger decrease of the condensate proportion; thus for an optically thick sample 
(where many reabsorptions take place) the number of condensed atoms would be reduced 
dramatically. However, these arguments are not rigorously valid in general. In ref. E3| it was 
reported that a fully quantum treatment shows that the reabsorption processes can, under 
certain conditions, help to increase the proportion of atoms in the BEC. This counterintuitive 
result was explained by means of an interference effect between different paths (which include 
reabsorptions) that lead to the same final states. This effect offers an interesting mechanism 
of continuous refilling and loading of a condensate. 

The treatment of ref. ]T^ | was based on an extremely simplified ID model in which all 
the ground states were taken into account, but in which there was a single excited level. In 
addition, the calculations of ref. jf2| did not consider the different Frank-Condon factors for 
the different possible decays. In this Letter, we analyze the BAR regime in a much more 
general model. We consider a 3D trap, and the possibility of more than one excited trap level. 
In addition, the Frank-Condon factors for the different decays are explicitely included in the 
calculations. We demonstrate that even for this much more complicated situation, under some 
conditions the reabsorption still helps to load atoms into the condensate. 

Let us consider a set of bosonic atoms with two internal levels |<?) and |e) confined in 
a dipole harmonic trap, which for simplicity is considered isotropic, with frequency u). We 
denote as N m the population of the m-th level of the \g) trap, where m = (m x , Tn v , m z ), and 
assume that the ground state trap verifies the BAR conditions. We analyze the situation in 
which a single atom in some state |e, I) (I = (l x , l y , l z )) decays via spontaneous emission into 
the \g) trap, producing a photon which can eventually be reabsorbed by another \g) atom 
(in particular by a condensed one), which in turn can later decay again, and so on. At some 
finite time, the scattered photon is no more reabsorbed and leaves the system. It is our aim 
to investigate how the atom distribution in the \g) trap changes during this process. 

For simplicity we do not consider the collisional mean-field effects. For typical s-wave 
scattering lengths this implies that in principle our results could be obtained for few atoms 
only (since we consider relatively small traps); one should point out , however, that although 
our calculations are limited to small traps, the BAR effect should be present also for larger 
ones, and in the presence of collisions. Although the main purpose of this Letter is to discuss 
a fundamental quantum effect and to present the methodology needed to study it, we must 
also point out that recent experiments |H| have shown that the s-wave scattering length can 
be modified by using Feshbach resonances. In particular, the regime of quasi-ideal gas can be 
experimentally achieved, in which the mean-field energy can be considered smaller than the 
trap energy. In such regime, the collisions introduce just a thermalization mechanism pi], 
and two and three-body collisional losses are almost absent |1| . For such modified scattering 
length our results are valid for much larger N . 
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Starting from the Hamiltonian which describes the bosons interacting with the quantized 
electromagnetic field, and using standard techniques, one can derive the master equation (ME) 
for the reduced density operator for the atomic degrees of freedom ]T^] , which assuming fi = 1 
becomes 

(1) 



where 



p = -iH e ffp + ipHl ff + Jp 



2=0 fc=0 

is a non hcrmitian effective Hamiltonian, and 



(3) 



1,1' m,m' 



is the jump operator. Here, ei (g m ) is the annihilation operator for atoms in the l-th (ru- 
th) excited (ground) level, u>f (uj^J is the energy corresponding to such state, and T is the 
spontaneous emission rate from the excited level. The complex coefficients ai mm >i> (whose 
explicit form can be found in @) are related to the Frank-Condon factors for the different 
excited-ground transitions. We denote in the following rji m the Franck-Condon factor for the 
transition between the Z-th excited level and the m-th ground level. The real part of the 
cuimm'V coefficient is denoted o§ mm , v . We consider an initial situation in which N m atoms 
occupy the m-th state of the ground-state trap, and a single excited atom is placed (with some 
probability given by a thermal distribution) in a state j of its corresponding trap. We denote 
the initial state as \ipo), and therefore the initial density matrix is defined as pjy = \^po)(ipo\. 
After a photon is released without further reabsorptions, we obtain the formal solution pn+i 
of the ME (|]), and calculate the probability to obtain a particular final state \ifjf) with N' m 
atoms in the m-th \g) trap level (and no atom in the excited state trap). After expanding 
into powers of the small parameter Nq 1 ^ 2 , and neglecting terms of order 0(Nq 3 ^ 2 ), such 
probability takes the form 



(ij>f\p N+ t\ipf) = — 

iVr 



o Jo 



dt I — 

4ir 
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with A (t) = exp[-iH^f f t] and A 1 = -i /* dr exp[-iH^ f (t - t)]H^\ exp[-iH^ f r] being 

1/2 

the terms of exp[— iH e fft] of order and 1 in 1/N Q ' , where 
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and 



H^ff = X! { a lrnWgl,gm + a l0m l> 5m50 }4 , 

1,1' m^O 
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are respectively the terms of H e ff of zeroth and first order in 1/N . Therefore, beginning 
from an initial relative number of condensate particles n — Nq/N, up to order 0(1/Nq), the 



4 



EUROPHYSICS LETTERS 



new relative number of condensed atoms after the last decay becomes 

,_ ^0 + 1 + ^0+2-^0) 
71 N + 1 



(7) 



where Pn +2 {Pn ) is the probability to have after the process Nq = No + 2 (No). Let \ip^' sS> ) 
be the final state with = N Q + 2, N' s = N s - 1, and A^ ,2 = ^Vj- The probability to decay 
into such a state is given by: 



N +2 



I / (2,s) 1 1 , (2,s)\ 



N 



dt 
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The corresponding process is schematically represented in fig. [l](a). Let us also define 1^'^) 
as the final state with Nq = N , N' s = N s + 1 and N'^q 2 = Nj . The probability to decay into 
this state takes the form: 



(0,a)\ 



N 



dt 



Jo 



n 
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where -R denotes the real part. The first and second terms in eq. (g) are respectively depicted 
in fig. 0(b) and (c), whereas the last term corresponds to the intereference between the en- 
ergetically equivalent paths which are shown in that figures. From eqs. (||),@ one obtains 
Pn +2 — J2 S < Pno+2 >' an< ^ Pn = S s < Pn > > wncrc <> denotes the average of the 
possible initial ground-state populations following the corresponding Bose-Einstein distribu- 
tion. We have numerically calculated the probabilities Pn and Pn q +2- From eqs. it 
becomes clear that such calculation requires to take into account all the possible paths con- 
necting a particular initial and final state which eventually involves an arbitrary number of 
emission-reabsorption cycles. This fact by itself makes the calculation extremely demanding. 
In addition, several important technical difficulties appear. First, we must note that in order 
to calculate the coefficients A$ and A\, one has to evaluate the exponential exp(— iMt), where 
M — v aiooi'ejeir is not an Hermitian matrix. Therefore, it becomes neccesary to introduce 

a biorthogonal set of operators /jj(fc) = J2i v ! i (^) e l > an( i /i(^) = Yli^l (^)e;, where vf'(k) 
and v[(k) are the right and the left (complex conjugate) eigenvectors of M. Therefore, such 
eigenvectors, and eigenvalues have to be calculated. Second, we note that the a coefficients 
constitute a 12 dimensional tensor, and that in order to evaluate every component, it is necce- 
sary to calculate the Franck-Condon factors for each possible excited-ground transition, and 
to perform the non-trivial task of calculating the Cauchy principal part integral appearing in 
the imaginary part of the a coefficients p2[ . Such strong technical difficulties make eventu- 
ally impossible to solve in 3D in a reasonable computational time systems with more than 4 
excited-state shells (20 levels), and therefore we have been constrained by such limit. 

We have analyzed the case of different number of energy shells of the e) and \g) trap, 
for different temperatures and different number of atoms. For each case, we have calculated 
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Fig. 1 - Schematic representation of the different processes appearing in eqs. (|8|) and (H). Straight 
lines denote atomic transitions, whereas waved lines represent the photons produced in the decays. 
Note the possibility of reabsorption in figs, (a) and (b). 



the probabilites (g|),(|9|), in order to evaluate the new relative number of condensed atoms 
n' . For the particular case of a single level in the excited state trap, results similar to those 
obtained in rcf. |Q were obtained. Figure ^ shows for different temperatures of the excited 
and ground trap, the case of 4 shells (20 levels) in the |e) trap, 10 shells (220 levels) in the 
\g) trap, N = 7 x 10 4 atoms, and a Lamb-Dicke parameter r/ 2 = uj r /uj = 2, where ui r is the 
recoil energy of the scattered photon. One can observe that for certain temperatures T e and 
T g of the |e) and \g) traps, the change in the relative number of condensed particles, n' — n 
is maximally positive, i.e. at these maxima the reabsorption effects help in the most efficient 
way to load the condensate. Away from such maxima, n! — n is positive for low temperature 
T e of the excited state trap and not very low temperatures T g of the \g). For very low 
T g , however, the positive processes described by P/vo+2 tend to vanish, since the number of 
non-condensed atoms which eventually could be repumped to the condensate becomes very 
small. On the top of fig. 0, such regions are those enclosed by the contours. Out of these 
regions, n' — n becomes negative, i.e. the reabsorption tends to decrease the condensation 
relative number. In all our calculations the BAR expansion has been proved to be valid, by 
checking the condition Pn q , Pn +2 <§C 1. We have in general observed that the BAR can fail 
not only for large T g , as expected from ref. |l2| ], but eventually also, for a given total number 
of atoms N, for large T e at the n' — n maxima. The latter is due to the large values of the 
imaginary part of the a coefficients at those peaks under such conditions, which invalid the 
expansion performed to obtain eqs. (^) and @. In the presented example, the condition above 
is fulfilled except for the maxima in the region ksTe > w and fcgT 3 > 50w. As discussed in 
ref. fl|], the phenomenon behind the positive effects of the reabsorption [pi — n > 0) cannot 
be explained by using (classical) rate equations, since it is given by the interference between 
the different paths which lead to the same final state. In particular, it becomes decisive that 
the interference term in eq. (||) is always destructive, since the process of fig. 2(b) includes 
an additional absorption-emission cycle (which gives a minus sign in the amplitude, like for 
a 27r-laser pulse). The consequence of this interference is that Pn decreases, which favours 
that the excited atom goes to the state \Nq + 1), and therefore it contributes to an increase 
of the proportion of condensed atoms. 

We have numerically simulated the pumping of atoms into a BEC in the \g) trap, under 
the BAR conditions. In our calculations, we have assumed that the collisions act on a shorter 
time scale than the spontaneous optical pumping. In that case, collisions will provide a fast 
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Fig. 2 - Change of the condensate fraction after a decay process as a function of the temperature 
T e of the |e)-trap and T g of the |g)-trap, the latter in exponential scale. The case of a total number 
of atoms N — 7 x 10 4 , 10 energy shells in the |g)-trap, 4 energy shells in the |e)-trap and r/ 2 — 2 is 
considered. On the top of the figure, the regions enclosed by the contours denote those temperatures 
for which the condensate fraction increases during the decay process 



thcrmalization mechanism of the \g) trap, between the two pump acts. After each pumping 
process we calculate the new condensate fraction, and accordingly the new T g , which we 
employ to evaluate the next pumping, and so on. In fig. |^ we analyze the same trap as in fig. 
||, for ksT e = u, and an initial number of trapped atoms N = 5 x 10 4 , for initial No/N = 0.9, 
0.98 and 0.99. One can observe that the positive effects of the reabsorption in the BAR regime 
allow to increase the condensate fraction during the pumping process, up to almost complete 
condensation. However, as pointed out previously, for very low T g the positive effects of BAR 
vanish. This explains the fact that the Nq/N is slightly lower than 1, and also that for the 
case of initial Nq/N = 0.99 the condensate fraction initially decreases. 

In this Letter we have extended the results presented in ref. Jl^] for a much more general 
3D situation, in which the excited state atoms can occupy more than one trapped state. We 
have shown that under the appropriate conditions, the BAR expansion is still valid for this 
more general case, and that the reabsorptions can play a positive role in the loading of the 
condensate. Such effect is a consequence of the always destructive interference between the 
processes which tend to lower the condensate fraction. The generalized treatment shows that 
the BAR condition presented in jl2| is not enough to guarantee the BAR expansion, since for 
more than one excited-trap level, the temperature of such trap is also important, and must 
be kept sufficiently low. Although for complexity reasons we have not analysed the situation 
in which the trap levels are distorted by the mean field provided by the atom-atom colli- 
sions, we must stress that similar analysis could be applied also if the atom-atom collisions 
modify the levels, if instead of the bare trap levels, self-consistent levels were considered, 
as in ref. |]r7[ . Therefore, under the BAR conditions, i.e. large condensation, the reabsorp- 
tion processes favour the optical pumping of atoms via spontaneous emission into a BEC. 
Such quantum effect is of fundamental interest, but can also be important experimentally to 
overcome condensate losses, and to achieve of a continuously-loaded atom laser. 
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Fig. 3 - Mean condensate fraction as a function of the number of loaded atoms, for the case of 10 
energy shells in the |<7)-trap, 4 energy shells in the |e)-trap and rj 2 = 2, and an initial total number 
of atoms N = 5 x 10 4 . From top to bottom, the curves represent, respectively, the case of an initial 
number of condensed particles No = 49632, 48828 and 45211. 
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